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Variabilidad estacional e interanual de la temperatura 
superficial en los mares interamericanos: patrones de 
sesgo en el sistema regional de modelacion del océano 


Abstract 

The Inter-Americas Seas (IAS), involving the Gulf of Mexico, the Caribbean and a 
section of the eastern tropical Pacific Ocean bordering Central America, Colombia and 
Ecuador, exhibits very active ocean-land-atmosphere interactions that impact socio- 
economic activities within and beyond the region, and that are still not well understood 
or represented in state-of-the-art models. On seasonal-to-interannual timescales, the 
main source of variability of this geographical area is related to interactions between 
Licencia Creative Commons the Pacific and the Atlantic oceans, involving anomalous sea-surface temperature (SST) 
DOT OUIEE §— oatterns like El Niho-Southern Oscillation (ENSO), and regional features in the Caribbean 
linked to the bi-modal seasonality of the Caribbean Low-Level Jet. This study investigates 
&) ‘O) © seasonal-to-interannual IAS surface-temperature anomalies in observations, and their 
representation in an eddy-permitting, 1/9° (or 0.11°; approximately 10 km) resolution 
simulation using the Regional Ocean Modeling System (ROMS), interannually-forced 
by the Climate Forecast System Reanalysis for the study period 1999-2008. Here, rather 
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Edited by: than analyzing model biases locally (i.e., gridbox-by-gridbox), a non-local SST pattern- 
René Parra based diagnostic was conducted via a principal component analysis. Regarding spatial 
Recibido / resolution, the results of this study suggest that even at eddy-permitting resolutions 
Received: of 0.11° several biases remain, which contrast with the general consensus in the 


a are literature. The approach permitted the identification of magnitude, variance and spatial 

Aceptado / systematic errors in SST patterns related to the Western Hemisphere Warm Pool, ENSO, 

fe aon the Inter-Americas Seas Dipole, and several other variability modes. These biases are 

mainly related to errors in surface heat fluxes, misrepresentation of air-sea interactions 

Bor ae eae! impacting surface latent cooling in the Caribbean, and too strong sub-surface thermal 
15/12/2021 stratification, mostly off the coast of Ecuador and northern Peru. 
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Los mares interamericanos (IAS), que se extienden desde el Golfo de México, el 
Caribe y la parte tropico-oriental del Océano Pacifico que limita con América Central, 
Colombia y Ecuador; presentan interacciones muy activas entre el océano, el suelo 
y la atmosfera, que influyen en las actividades socioeconémicas dentro y fuera de la 
region. Estas interacciones atin no se comprenden en su totalidad y no se representan 
adecuadamente, inclusive en modelos de ultima generacion. En escalas de tiempo 
estacionales e interanuales, la principal fuente de variabilidad en esta zona geografica, 
se relaciona con las interacciones entre los océanos Pacifico y Atlantico, incluyendo 
patrones andmalos de la temperatura en la superficie del mar (SST), como El 
Nifto-Oscilacién del Sur (ENSO), y otras caracteristicas regionales relacionadas con la 
estacionalidad bimodal del chorro de bajo nivel del Caribe. Este estudio investiga las 
anomalias estacionales e interanuales de los registros de la temperatura en superficie 
en los mares interamericanos, y de los resultados simulados con una resolucién 
a una escala de remolino, de 1/9° (0,11°, aproximadamente 10 km), generados 
por el Sistema Regional de Modelado del Océano (ROMS), y en base del Reanalisis 
del Sistema de Pronostico del Clima para el periodo 1999-2008. En lugar de analizar 
os sesgos de los resultados del modelo de manera local (cuadricula por cuadricula), 
se realizo un diagndstico no local, en base de la identificacidn de patrones de la 
SST, mediante un analisis de componentes principales. Los resultados indican que 
incluso con resoluciones a escala de remolino, y a diferencia del consenso de la 
iteratura, persisten varios sesgos. El enfoque permitid cuantificar la magnitud, la 
varianza y los errores sistematicos de los patrones espaciales de la SST, relacionados 
con la piscina calida del hemisferio occidental, ENSO, el Dipolo de los mares 
interamericanos y otros componentes de variabilidad. Estos sesgos se relacionan 
principalmente con errores la cuantificacién de los flujos de calor en la superficie, 
una inadecuada representacién de las interacciones entre el mar y la atmésfera, 
componentes que afectan el enfriamiento latente en superficie en la zona del Caribe 
y genera una estratificacién térmica subterranea excesiva, principalmente frente a la 
costa de Ecuador y el norte de Peru. 


Palabras clave: mares interamericanos, Sistema Regional de Modelizacién del Océano 
(ROMS), evaluacion de modelos en base de patrones. 


INTRODUCTION 


Sea-surface temperatures (SST) are key to understanding physical mechanisms 
regulating regional and global climate at multiple timescales. In the Inter-Americas 
Seas (IAS), normally referred to as “Intra-Americas Seas,’ a region involving the Eastern 
Tropical Pacific, the Gulf of Mexico and the Caribbean, SST patterns are subject to 
large spatial and temporal biases in state-of-the-art modeling systems, especially in 
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imitations in the availability of this kind of data 


is through the use of the numerical ocean models. Although several model simulations 
exist, few evaluations exist of how well these products represent key oceanic features in 


this complex region. 


Among the numerical models used in the study of ocean circulation processes at 
regional level, the Regional Ocean Modeling System (ROMS) stands out [12-13]. This 
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model solves the hydrostatic, incomprehensible fluid primitive equations using "terrain- 
following" curvilinear vertical coordinates [14]. ROMS has been used by the Colombian 
Navy since 2012 following the operational forecast configuration of the Latin American 
Observatory for Climate Events [15], but this is the first report of the formal evaluation 
conducted on the model. 


The main objective of this study is to identify how well oceanic conditions, specifically 
sea-surface temperature patterns, are represented in a realistically-forced interannual 
simulation in ROMS in the IAS. To accomplish this goal, a pattern-based evaluation 
approach is reported, focusing on the analysis of observed and simulated key modes of 
SST variability along the year. In the following section the region of study is described. 
Data and methods are introduced in Section 3, followed by the main results in Section 
4. The overall discussion of the results is presented in Section 5, and the concluding 
remarks are summarized in Section 6. 


STUDY AREA: THE INTER-AMERICAS SEAS 


The Inter-Americas Seas is part of the WHWP, defined as a broad area including the Gulf 
of Mexico, southern United States, Mexico, the Caribbean Sea with its islands, Central 
America, northern South America, and the ocean off the west coast of Central America 
and Colombia (Fig. 1). 
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FIGURE 1. Geographical area of study in the Inter-Americas Seas. Colors represent depth in meters below mean 
sea level. 
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According to some authors, the WHWP is related to the genesis of the warm and cold 
phases of ENSO [16-17]. Between the boreal summer and autumn, the WHWP exhibits 
a peak in SST [18-19] that dominates most of the eastern tropical Pacific region [20- 
21], crossing Central America to encompass areas of the Caribbean and the Central 
Atlantic Ocean. 


During winter (December to February, DJF), the Inter-Tropical Convergence Zone (ITCZ) is 
at its most southern position, and the northeastern trade winds intensify in the region [22]. 
The difference in atmospheric pressure between the Gulf of Mexico and the Caribbean, 
as well as in the Tropical Pacific Ocean (TPO), is related to the presence of the Caribbean 
Low-Level Jet (CLU, [23-24]) and the bi-annual incursion of NASH in the region of interest 
25], with intense winds that cross the mountainous area of Central America and interact 
with other low-level jets in the gulfs of Tehuantepec, Papagayo and Panama [26-27]. The 
Papagayo Jet crosses the depression of Lake Nicaragua and rises on the Pacific coast, 
about 70 km north of the Gulf of Papagayo, near San Juan del Sur, Nicaragua. The Panama 
Jet extends from the coastal area of the Gulf of Panama, centered around 79°W, to the 
Ecuadorian maritime territory. As a result of the interaction of these jets with the upper 
ocean layer, there is a coastal upwelling (Ekman positive pumping) that generates a zonal 
gradient in the SST in the Panama Bight region, as well as in the maritime territory of 

exico and Nicaragua in the Pacific Ocean [23-24, 28]. In general, the SST isotherms over 
he Caribbean and the eastern tropical Pacific are mostly zonally distributed. 


DATA AND METHODS 
Interannual simulation using the Regional Ocean Modeling System (ROMS) 


In the present simulation, the lateral boundary conditions used were taken from the 
Simple Ocean Data Assimilation (SODA) dataset [29]. The surface forcing of the model 
was carried out with data from the National Center for Environmental Prediction - 
Climate Forecast System Reanalysis (NCEP-CFSR), which provides global atmospheric 
conditions from 1979 to 2010 [30]. These data have a spatial resolution of approximately 
38 km (7382), with global coverage and 64 pressure levels. 


The Extended Reconstructed Sea-Surface Temperature dataset version 5 (ERSSTV5) 
provides sea-surface temperature fields at 2 degrees and monthly resolutions for the 
entire planet, from January 1854 to present [31]. 


ROMS is a high-resolution model that uses curvilinear coordinates that follow the 
bathymetry and solves the primitive equations in a rotating reference system (the 
Earth), based on the Boussinesq approximation and considering hydrostatic [32]. |n this 
study, version 3.0 of ROMS AGRIF was adopted, with which an interannual experiment 
was carried out through a 12-year period. This model has been implemented previously 
for the studies in the Southeast Pacific, analyzing the mean circulation, the seasonal 
cycle, and the mesoscale dynamics of the Peru Current System (PCS; [13]), exhibiting 
a good representation of the predominant ocean features in this sector. Nonetheless, 
a pattern-based evaluation of the representation by ROMS of SST patterns in the IAS 
region has not been conducted until now. 
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Model configuration 

In the implementation of the model, a computational mesh was configured that 
avanees. frames a section of the IAS included in the region defined by 100°W to 56°W, and 
mene 10°S to 25°S. The bathymetric information has a resolution of two (2) arc-minutes and 

was obtained from the global ETOPO2 dataset, which uses observations of depth 

sounding and satellite altimetry [33]. 


Although ROMS has terrain-following coordinates that allow for a more realistic 
representation of the bathymetry, a smoothing of the terrain slopes is used to avoid 
errors in the pressure gradient [34]. This smoothing is controlled in the model via the 
rtarget parameter [35], which, after some sensitivity experiments, was assigned a value 
of 0.15. 


The resolution of the simulation is 1/9° (or 0.11°; approximately 10 km) forthe construction 
of the horizontal grid, and 32 levels of depth for the vertical sigrna coordinate. For the 
distribution of the vertical layers, values of surface stretching and bottom stretching, 
respectively, 8, = 6 and 9, = 0 were used, corresponding, respectively, to the coordinate 
parameters for the surface and bottom levels. 


A total of 12 years were simulated, from January 1997 to December 2008. By the end of 
the second simulated year, the volume-averaged kinetic energy reached an asymptotic 
value. Hence, the first two years were discarded due to the model's spin-up period, and 
the remaining 1999-2008 period was used in the analysis. Two simulation years have also 
been reported by previous studies (e.g., [36, 13]) as a typical time for model stabilization. 


Data processing and analysis 


The mean seasonal cycle was first subtracted from the sea-surface temperature fields 
for each non-concurrent season of interest in the period 1999-2008: DJF, MAM, JJA 
and SON (respectively periods of every three months, December to February, March 
to May, June to August and September to November). Then a Principal Component 
Analysis (PCA) was conducted on the sea-surface temperature anomalies, SSTa, for 
each season using the International Research Institute for Climate and Society's 
Climate Predictability Tool (CPTv16; [37]), which computes the Empirical Orthogonal 
Functions via a singular-value decomposition. The first three Empirical Orthogonal 
Functions (EOFs) were used for the analysis, typically explaining 73-90% of the total 
observed variance, depending on the season. 


Since the sample size for the PCA is relatively short, a comparison was conducted 
between the observed EOFs obtained when using the same 10-year period of the 
simulation (1999-2008), and a 20-year (1994-2014) and a 30-year (1984-2014) period, 
which included strong El Nifo and La Nifa events. No major differences were obtained 
(not shown) in the observed patterns or the explained variance associated with each 
EOF in any of the seasons, and hence the selected 10-year period was considered 
enough for the analysis. 


Because the spatial horizontal resolution of the simulation does not match the 
horizontal resolution of the ERSSTv5 dataset, the latter was interpolated from ~0.11° to 


6 _ DOI: https://doi.org/10.18272/aci.v13i2.2046 


Articulo/Article 
Seccion/Section B 


Vol. 13, nro. 2 
ID: 2046 


avances 
en ciencias e 
Ingenieria 


7 


Seasonal-to-interannual variability of sea-surface temperatures in the inter-americas seas: 
pattern-dependent biases in the regional ocean modeling system 
Caicedo-Laurido / Mufioz / Chourio / Tobar / Latandret (2021) 


2°. Nonetheless, to avoid losing spatial features when computing the PCA, the analysis 
was also conducted using the original model resolution. 


Spatial patterns in the simulation (with and without interpolation) were compared 
against the observed EOFs for each season, and whenever major differences arose, model 
biases were related to concrete physical modes of variability and misrepresentation of 
physical mechanisms in the region. 


SEASONAL AND INTERANNUAL VARIABILITY OF SEA-SURFACE 
TEMPERATURE ANOMALIES IN THE IAS 


This section presents the main results of the pattern-based diagnostic analysis, in 
conjunction with the identification of biases in ROMS based on an extensive literature 
review for discussion. 


DJF season 


fog b 

~~ EOF3 
15% RB \ 

00°wW =90°W so°w 70°W 60°W 100W 8=(90°W so°w 70°W 60°W 


09 06 -03 00 030 060 0.90 


FIGURE 2. Spatial distribution of the first three modes (rows) of SSTa variability for the DJF season for ERSSTv5 
(a, d and g), ROMS (b, e and h) and ROMS interpolated at a resolution of 2° (c, fand i). Each panel displays the 
explained variance associated with each EOF. 
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The SSTa spatial pattern identified via the first Empirical Orthogonal Function (EOF1; 


variance exhibited by the simulated pattern is ~11% less 


han the observed one. 


Fig. 2a) sketches the seasonal imprint of the WHWP covering most of the analysis 
domain, consistent with the so-called Meridional-Wind Mode [38]. The EOF1 is related 
o northerly wind anomalies linked to below-normal SST in the Caribbean, and at the 
same time southeasterly wind anomalies over the Gulf of Mexico, coupled with above- 
normal SST there (Figs. 2a-c). This pattern is also associated with frequency variations of 
he high-latitude cold air-mass intrusion [39-40]. In this season, ROMS largely captures 
he main spatial features of the observed pattern (Fig. 2b, c), although the explained 


The EOF2 (Fig. 2d) exhibits a meridional dipolar structure in the Caribbean region, known 


as the IAS Dipole, that is mostly driven by changes in 


he air-sea heat fluxes [17]. 


A 


weak structure appears mainly in the maritime zone of Ecuador in the observed EOF2, 
contrasting with the rest of the eastern Pacific. Contrary to the performance exhibi 


by ROMS on reproducing the main features of the WH 
he model erroneously interchanges EOF2 and EOF3 (s 
in Fig. 2), assigning slightly more (~4%) explained va 


iance than it should to 


hird observed EOF. Furthermore, the simulated EOF3 (cfr. Fig. 2d, h-i) shows bia 


both in the magnitudes and geographical extension o 


west of Cuba. 


in a future study. 


surface thermal stratification. 


o the northeastern sector of the domain, showing just a weak signal over and to 


ed 


WP seasonal pattern (EOF 1), 
ee middle and bottom rows 


he 
ses 


f the IAS Dipole, confining it 


he 


n addition, ROMS overdoes the observed SSTa pattern over the equatorial Pacific 
cfr. Fig. 2d, h-i), suggesting a warm bias related to a misrepresentation of the 
weak sub-surface thermal stratification (“thermostad”) observed in the region [41] 
and reported in certain ROMS configurations [42]. The higher spatial resolution of 
ROMS allows resolving the effects on SSTa of the Tehuantepec, Papagayo and Panama 
Low-Level Jets in the Pacific Ocean region in the observed EOF2, an influence that 
is also observed in EOF3, especially in the case of the Panama Jet. Nonetheless, their 
impact on SST fields cannot be formally evaluated in this work, as ERSSTv5 does not 
resolve them adequately. This topic deserves further analysis and will be discussed 


The Caribbean section of the observed EOF3 (Fig. 2g) is related to the seasonal 
intensification of the CLL, bringing cooler-than-average SSTa off the coast of northern 
South America (Fig. 2e, f), in part contributing to rainfall inhibition during this season. 
Except for the fact that ROMS confused EOF2 and EOF3, the magnitudes and spatial 
patterns are represented by the model, although the SSTa pattern off the coast of 
Ecuador is stronger than it should be, suggesting again model biases related to sub- 


The actual difference in the explained variance of the simulated EOF2 and EOF3 is 
~1% (cfr. Fig. 2e, h), explaining the inversion in the ordering of the simulated EOFs. 
Nonetheless, the same difference in the ERSSTV5 dataset is ~5% (Fig. 2d, g). Even when 
these differences can be considered small, they are important when analyzing the 
modes of variability that most contribute to explain the SSTa variance and indicate 
biases in the physical processes represented in the model. 
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FIGURE 3. Spatial distribution of the first three modes (rows) of SSTa variability for the MAM season for ERSSTv5 
(a, d and g), ROMS (b, e and h) and ROMS interpolated at a resolution of 2° (c, f and i). Each panel displays the 
explained variance associated with each EOF. 


In this season, the first observed spatial pattern explaining most of the SSTa variance 
(EOF1, Fig. 3a) is consistent with the WHWP SST distribution in both the Pacific and 
the Caribbean. Although ROMS assigns a lower explained variance (~7% less than the 
observed one; cfr. Fig. 3a-c) to this mode of variability, it tends to reproduce the large- 
scale features observed, related to a homogeneous pattern of SST between 27°C and 
28°C present in the Caribbean Sea during the season [43]. 


EOF2 (Fig. 3d) shows a spatial distribution associated with an IAS [17, 38], which is 
characterized by late-winter and early-spring positive SSTa in the Caribbean Sea and 
negative SSTa in the Gulf of Mexico. This dipole pattern is not really present in the ROMS 
simulation (Fig. 3e-f), suggesting biases related to the representation of surface heat fluxes 
over the Caribbean and the Gulf of Mexico. Nonetheless, the simulation assigns to this mode 
a value of the explained variance that is similar to the one observed (a difference of ~2%). 
This relatively high explained variance in the simulation is related to the model assigning a 
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large statistical charge to the spatial pattern in the eastern tropical Pacific, related to ENSO, 
which otherwise exhibits a similar spatial distribution to the ERSSTV5 (Fig. 2d-f). 


The observed EOF3 (Fig. 3g) corresponds to the “Transition Mode” [38], exhibiting a 
dipolar pattern in the Caribbean but contrasting with the IAS Dipole (EOF2) because the 
signal is almost absent in the Gulf of Mexico, and does not extend into the Pacific, which 
shows a mostly homogeneous spatial distribution of SST. Compared to the observed 
EOF3, ROMS exhibits in the Pacific a zonally-distributed SSTs pattern of opposite sign 
along 5°N and also a noisy SSTa region off the coast of Ecuador and northern Peru; 
similarly, the simulation shows a dipolar structure too strong and extending farther 
west into the Gulf of Mexico, compared to the observed pattern (Fig. 3h, i). Furthermore, 
ROMS underestimates the explained variance for this mode of variability in ~6%. 


In summary, in the boreal spring ROMS tends to represent well the spatial features of 
the ENSO peak season in the Pacific, and the WHWP in the IAS, but there are important 
biases regarding the representation of the IAS Dipole. 
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FIGURE 4. Spatial distribution of the first three modes (rows) of SSTa variability for the JJA season for ERSSTv5 
(a, d and g), ROMS (b, e and h) and ROMS interpolated at a resolution of 2° (c, f andi). Each panel displays the 
explained variance associated with each EOF. 
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WP (EOF 1; Fig. 4a) pattern has expanded, covering a larger 


extension than in the previous season in both the Caribbean and the Pacific [18, 44- 


45]. The WHWP normally tran 


sitions into the AWP in June and shows a peak typically 


in early September [1]. This spatial pattern is better represented by ROMS (Fig. 4b, c) in 


the Caribbean than in the Paci 


fic. Although similar patterns appear in both ERSSTV5 and 


ROMS off the Pacific coast of Central America, in ROMS the pattern is noisier than the 
ne observed for the rest of the Pacific region (Fig. 4b, c). The difference in explained 


(e) 


variance between 


Most of the signal 


he observed and simulated EOF1 is ~12%. 


in the observed and simulated EOF2 (Fig. 4d-f) appears in the Pacific 


region, which exhibits a pattern associated with ENSO-like SSTa [46-48]. The simulation 


and observed exp 


~1%), indicating ROMS tends to reproduce adequately 


spatial distribution of SSTa variability during th 


The EOF3 (Fig. 4g-i) resembles the signature o 


n 


h 


Vv 
p 


in ERSSTv5 and with a slightly more equatorial 


ia a strengthening of the zonal SST gradient i 
romoting stronger upwelling and thus creating a positive feedback that reinforces 


ode [49], with its narrow zonal SSTa pattern 
orthern Peru. This mode, which ROMS reprod 


igher resolution of the simulation—is repor 


the SST contrast. A homogeneous SSTa patte 
simulated EOF3 for most of the Caribbean, with an opposite sign 


the so-ca 


uces with 


ed to be 


ained variances for this mode are basically the same (difference of 


he mechanisms related to the 


e mature season of ENSO. 


led eas 


tern Pacific SST Cooling 


extending off the coast of Ecuador and 


he same explained variance as 


elated 


nis observed in 


In summary, in the boreal summer the model has a fair represen 
spatial features of the modes analyzed, and the total explained variances the EOF2 and 
EOF3, except for the considerably less homogenous signal of the WHWP in the Pacific in 
EOF1, which consistently exhibits less variance in ROMS than in 
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FIGURE 5. Spatial distribution of the first three modes (rows) of SSTa variability for the SON season for ERSSTv5 


(a, d and g), ROMS (b, e and h) and ROMS interpolated at a resolution of 2° (c, f and i). Each panel displays the 
explained variance associated with each EOF. 


As with the other seasons, most of the observed explained variance (~47%) is associated 
with the WHWP (EOF 1, Fig. 5a), which exhibits very low SSTa variability in the Caribbean 
off the coast of Venezuela and Colombia during the SON season [18], and during 
November involves the AWP demise phase [1]. Although ROMS does capture most of 
the spatial features (Fig. 5b-c), especially after interpolating at 2°, it assigns an explained 
variance which is ~13% lower than the observed one. 


The second EOF (Fig. 5d-f) corresponds to the Interocean Mode [38, 1] and is related 
to southwesterly winds that contribute to the formation of the inter-basin dipolar 
configuration. The pattern also has a contribution related to the second annual weakening 
of the CLL, which favors an increase in positive SSTa [24] between ~12°N-17°N and 80°W- 
70°W. The Interocean Mode tends to be enhanced during ENSO events, which modulate 
the strength of the CLU, reinforcing the SSTa dipolar configuration [38, 50]. The main 
observed spatial feature of the SON EOF2 in the tropical Pacific region resembles an 
ENSO-type configuration, consistent with previous studies (e.g. [51]). ROMS captures the 
main spatial features of the inter-basin dipole, but the explained variance is ~4% higher 


DOI: https://doi.org/10.18272/aci.v13i2.2046 


Articulo/Article 
Seccién/Section B 


Vol. 13, nro. 2 
ID: 2046 


avances 
en ciencias e 
Ingenieria 


13 


Seasonal-to-interannual variability of sea-surface temperatures in the inter-americas seas: 
pattern-dependent biases in the regional ocean modeling system 
Caicedo-Laurido / Mufioz / Chourio / Tobar / Latandret (2021) 


han that of the observed EOF2. Part of this extra variance can be related to the fact that 
ROMS tends to represent high-resolution features, like the Papagayo Jet (southwest off 
he coast of Nicaragua), which are less notorious in the ERSSTv5 dataset (Fig. 5d-f). Other 
differences are related to a local pattern off the coast of Ecuador in the simulation, that is 
not really visible in ERSSTv5; and a mostly inverted signal, and noisier north of 20°N than 
he observed pattern, in the Caribbean (Fig. 5d-f). 


The Caribbean part of the EOF3, both in ERSSTv5 and ROMS (Fig. 5g-i), resembles the 
AS Dipole discussed in the analysis for the DJF season (section 4.1). EOF3 is considered 
here a precursor of the IAS Dipole, related to the so-called Meridional-Wind Mode [38]. 
Although the difference in explained variance between the ERSSTv5 and ROMS in this 
EOF is ~1%, the simulation tends to amplify both in magnitude and areal extension the 
observed SSTa configuration along the Pacific coast of Central America. 


nsummary, in the boreal fall the main model spatial-pattern biases are related to a noisy 
SSTa field north of 20°N (EOF2), and a too strong and too wide SSTa signal in the Pacific, 
north of the Equator and south of the Central American coast (EOF3). 


Main patterns and biases in ROMS 


The most common approach when evaluating model biases is to report a variety of 
statistical metrics to quantify the differences between model output and ERSSTV5. 
Although this approach is deemed very useful by the scientific community, a drawback 
is that it tends to lack a direct physical interpretation of the reported errors. For example, 
frequently used metrics like Pearson's correlation or Root Mean Square Error can be very 
useful to quantify how in-phase two variables are, or an average measure of how different 
he magnitudes of those variables are, respectively, but in and of themselves these metrics 
do not provide information about which physical processes are not being well represented 
in the model universe. 


Furthermore, the most common approach evaluates biases on a gridbox-by-gridbox 
or local) basis, thus making it difficult to assess if particular observed spatio-temporal 
patterns are being represented in the model. These patterns can in general be shifted, 
rotated, deformed, rescaled or actually not represented at all in the model world [52]. 


An alternative approach, used in this study, focuses on the evaluation of the model 
representation of the most important observed modes of variability of the field of 
interest—SST in this case—and that helps identify problems in the representation of key 
physical processes. In brief, a pattern-based evaluation approach permits understanding 
when the model reproduces the right processes for the right reasons and adds value for 
both scientific and operational purposes. 


Using a Principal Component Analysis to conduct the model pattern-based evaluation 
has the additional advantage of being able to use a Principal Component Regression 
(PCR) and other EOF-based methods (see [53]) as bias-correction approaches, also 
known as calibration or Model Output Statistics [54]. To make the case, consider PCR: 


1-28 (a, EOFx) + Bi (1) 
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Equation (1) builds a linear regression model for the value Y, of the variable of interest (SST 
in this case), for each gridbox i, where m counts for the maximum number of the model 
patterns--EOFx--considered (e.g., m=3 in this study), a, is the slope regression coefficient 
for the gridbox jand--EOFx--, and is the corresponding intercept value. In other words, a 
PCR is a multiple linear regression which uses the model EOFs as predictors rather than 
each gridbox’s value of the variable under consideration. 


Hence, analyzing how well the observed patterns are reproduced by a model also 
helps to better understand how much of the total bias is associated with each pattern, 
assessed, for example, via the explained variance reported by each EOF. This is the 
approach followed in this study, and the EOF analysis and associated variances were 
described in the previous section. 


The simulated patterns do not need to be perfect to be able to conduct a bias-correction 
through these EOF-based methods. If the most important observed modes of variability 
are simulated, even with biases that shift their location or rotate them, these methods 
are in principle able to extract useful information from the model. The Model Output 
Statistics approach is especially useful when interested in prediction (for an example 
considering the IAS, see [6]). Although this study does not directly deal with calibration 
and prediction, these topics underscore the importance of the analysis conducted here 
regarding better understanding of model biases in the main patterns of variability in 
ROMS, at least for SSTa in the IAS region. 


Overall, the main biases are summarized in Table 1. ROMS tends to correctly reproduce 
the main mode of variability (EOF1) in IAS, characterized by the WHWP signature, in 
all seasons but the boreal summer. Coincidentally, the second and third modes of 
variability tend to exhibit different errors in different seasons, but these tend to be less 
pronounced in the boreal summer. As indicated in the Introduction, most state-of-the- 
art global models tend to show their largest biases in boreal summer and fall [1-5]. 


TABLE 1. Main patterns and biases in ROMS. 


Season EOF Patterns Biases (w.r.t. ERSSTv5) 


Western Hemisphere Warm Pool 
1 WHWP; [18]) No major biases. 
Meridional-Wind Mode [38] 


ROMS erroneously interchanges 
EOF2 and EOF3. Biases in magnitude 
and extension of the IAS Dipole, 
suggesting errors in surface heat 
fluxes. SSTa warm biases off the coast 
of Ecuador related to too strong sub- 
surface thermal stratification. 


IAS Dipole [17] 
DJF Thermostad east of Galapagos [41] 


SSTa warm biases off the coast of 
Ecuador related to too strong sub- 
surface thermal stratification. Also 
biases along the Pacific coast of 
Central America. 


Caribbean Low-Level Jet 
(CLL; [23, 25]) 
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Season EOF Patterns Biases (w.r.t. ERSSTv5) 
1 WHWP [18] No major biases, 
Biases related to the representation 
yi IAS Dipole [17] of surface heat fluxes over the 
ENSO-like configuration [51] Caribbean and the Gulf of Mexico, 
MAM ENSO signal is too strong. 
SSTa biases off the coast of Ecuador 
5 and northern Peru. Dipolar SSTa 
Transition Mode [38] structure too strong and extending 
farther west into the Gulf of Mexico 
than in ERSSTV5. 
SSTa pattern represented by ROMS is 
| WHWP [18] noisier than the one observed in the 
Pacific region. 
He 2 ENSO-like SSTa [55]. No major biases. 
3 Eastern Pacific SST Cooling Mode No major biases. 
[46] 
Biases related to a SSTa local pattern 
1 WHWP [17] off the coast of Ecuador in the 
simulation. 
SSTa biases north of 20°N, potentially 
Interocean Mode [38] : : 
SON 2 ENSO-like configuration [18] related to misrepresentation of local 
surface latent cooling. 
ROMS tends to amplify both in 
fete scar magnitude and areal extension the 
: MalicionaWindiicac tee observed SSTa configuration along 
the Pacific coast of Central America. 
Generally speaking, ENSO-like and Pacific SSTa configurations—although sometimes 


too strong—tend to be better represented than some patterns in the Caribbean 


and Gulf of Mexico, with some key biases in the IAS Dipo 


e between boreal winter 


and spring, and SSTa patterns off the coast of Ecuador and northern Peru (winter), 


and alongshore the Pacific coast of Central America (fal 


| and winter). Consistent 


with previous studies [1], these errors are suggested here to be related to biases in 
surface heat fluxes, misrepresentation of the impact of westerly anomalies on surface 
latent cooling in regions of the Caribbean, and warm biases related to too strong 
sub-surface thermal stratification. Is it also possible that a component of the biases 
is related to the assimilation of initial and boundary conditions, as reported by some 
authors (e.g., [42)). 


Regarding spatial resolution, the results of this study suggest that even at eddy- 
permitting resolutions of 0.11° several biases remain, which contrast with the general 
consensus in the literature [1, 55-56]. It is also true that some spatial features cannot be 
adequately evaluated using ERSSTv5 due to the relatively low-resolution of the dataset, 
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such as the effects of the Tehuantepec, Papagayo and Panama low-level jets in the Pacific 
Ocean region. This is another reason why a process-based model evaluation approach 
is an important complement, as more traditional bias assessment between model and 
ERSSTv5 when the horizontal resolution is very different can be misleading, even after 
interpolation to a common horizontal resolution. 


CONCLUDING REMARKS 


This study investigated variance and spatial biases in sea-surface temperature patterns 
in the Inter-Americas Seas, represented by a 0.11° resolution, interannually-forced 
simulation using the Regional Ocean Modeling System. Rather than using a variety 
of statistical metrics to assess local (i.e, gridbox-by-gridbox) differences between the 
simulation and data based on observations (ERSSTv5), the study focused on analyzing if 
the sea-surface temperature patterns that explain most of the total observed variance 
are present in the simulation, and, if so, which are the main biases exhibited by each 
pattern. Furthermore, the study analyzed the most likely physical processes involved in 
those pattern-dependent biases. 


The large-scale spatial features and variances related to the Western Hemisphere Warm 
Pool and ENSO-like patterns tend to be well represented in the model in most but not 
all seasons. Overall, spatial biases are more frequently present in the Caribbean and 
Gulf of Mexico than in the eastern Pacific region under study, with magnitude, variance 
and areal extension biases related to the Inter-Americas Seas Dipole (boreal winter 
and spring), and the Transition (spring), Interocean and Meridional-Wind (fall) modes. 
The analysis suggests that model biases are mainly related to errors in surface heat 
fluxes, misrepresentation of air-sea interactions impacting surface latent cooling in the 
Caribbean, and too strong sub-surface thermal stratification, typically off the coast of 
Ecuador and northern Peru. 


The fact that several pattern-dependent biases exist in the eddy-permitting, relatively 
high-resolution simulation conducted suggests that important errors remain even in 
configurations running at high horizontal resolution and that additional work is needed 
to improve model representation of physical processes and adequate assimilation of 
initial and boundary conditions. 
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